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We provide an interpretation of the recent ATLAS diboson excess in terms of a class of supersym¬ 
metric models in which the scale of supersymmetry (SUSY) breaking is in the few TeV range. The 
particle responsible for the excess is the scalar superpartner of the Goldstone fermion associated 
with SUSY breaking, the sgoldstino. This scalar couples strongly to the Standard Model vector 
bosons and weakly to the fermions, with all coupling strengths determined by ratios of soft SUSY 
breaking parameters over the SUSY breaking scale. Explaining the ATLAS excess selects particular 
relations and ranges for the gaugino masses, while imposing no constraints on the other superpartner 
masses. Moreover, this signal hypothesis predicts a rate in the Z 7 final state that is expected to be 
observable at the LHC Run II already with a few fb“^ of integrated luminosity. 


I. INTRODUCTION 

The ATLAS Collaboration recently published a search 
for resonances in the boson tagged di-jet mass distribu¬ 
tion, featuring an excess of events around 2 TeV [1]. De¬ 
spite the fact that the statistical significance of the ex¬ 
cess (up to 3.4(7 locally and 2.5cr globally) is limited, the 
appearance of other excesses, though less significant, in 
similar final states and in the same mass region, such as 
in the CMS search in ref. [2], motivates some theoretical 
effort to understand the possible origin of these fluctu¬ 
ations. Several papers have already appeared, aimed at 
explaining the excess in terms of different new physics 
models [3-36] and scrutinizing the ATLAS analysis [37] . 

In the context of supersymmetry (SUSY), it is not 
straightforward to find an explanation of this excess. For 
instance, with the usual particle content of the mini¬ 
mal SUSY extension of the Standard Model (MSSM), 
there is no candidate particle that could give rise to such 
a signal. However, in the case where SUSY is broken 
at a low scale^, additional degrees of freedom, related 
to the spontaneous breaking of SUSY, are present and 
can become phenomenologically relevant. In particular, 
the Goldstone fermion of SUSY breaking, the goldstino, 
and, when SUSY is linearly realized, its scalar superpart¬ 
ner, the sgoldstino, can couple strongly to some of the 
SM particles. The interaction strengths of the goldstino 
and sgoldstino are determined by ratios of the usual soft 
SUSY breaking parameters of the MSSM over the super- 
symmetry breaking scale. 

In this paper we provide an interpretation of the AT¬ 
LAS diboson excess in terms of a class of SUSY models 
where the SUSY breaking scale is in the few TeV range, 
with a 2 TeV sgoldstino scalar being responsible for the 


^ The current experimental lower bound on the SUSY breaking 
scale is at or below 1 TeV [38] (the exact value depends on the 
superpartner spectrum). 


excess. For different discussions concerning sgoldstino 
physics, see, for example, refs. [39-49]. The sgoldstino 
couples mostly to the SM vector bosons, with interaction 
strengths determined by ratios of the gaugino masses over 
the SUSY breaking scale, whereas its couplings to the SM 
fermions are generically suppressed. In what follows, we 
study the compatibility of this signal hypothesis with the 
excess, identify the relevant region of the parameter space 
(in terms of the gaugino masses and the SUSY breaking 
scale) and discuss the relations to other searches in cor¬ 
related channels, such as 77 and Z^. 

The paper is organized as follows. In Section II we 
provide the sgoldstino couplings to the SM and in Sec¬ 
tion III we extract the relevant values for the gaugino 
masses and SUSY breaking scale that allow us to explain 
the ATLAS excess. We describe the constraints from, 
and implications for, other searches in Section IV and 
conclude in Section V. 

II. THE SGOLDSTINO MODEL 

If SUSY is realized in Nature, since the SM particles 
are not mass-degenerate with their superpartners, it must 
be in a broken phase at low energies. A general conse¬ 
quence of the spontaneous breaking of (global) SUSY is 
the existence of a Goldstone fermion, the goldstino. We 
will assume that the goldstino resides in a gauge singlet 
chiral superheld, with SUSY linearly realized, 

Y = X -k V2eG + 9^Fx (1) 

where the auxiliary field acquires a vacuum expectation 
value (VEV), {Fx)=f, that gives the dominant contri¬ 
bution to supersymmetry breaking. The focus of this 
paper will be on the complex scalar superpartner of the 
goldstino, the sgoldstino x in eq. (1). In contrast to the 
goldstino, the sgoldstino is not protected by the Gold¬ 
stone theorem and therefore it will generically acquire 
a mass, with a value that is model-dependent. Also, in 
general, the masses of the GP-even and GP-odd scalars 
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do not need to be equal. Here we assume them to be 
equal and fix them to be 2 TeV. 

One way to take into account the interactions of the 
goldstino and sgoldstino is to simply promote all the 
usual MSSM soft terms to SUSY operators involving the 
goldstino superfield in eq. (1). For instance, the gaugino 
masses m^, where i = 1, 2 and 3 corresponds to the bino, 
wino and gluino masses, respectively, are promoted to 
the following SUSY operators. 


^ iy / , ( 2 ) 

where for i=l,2 and 3, corresponds to the gauge 
field-strength superfield for U{1)y, SU{2)l and SU{3)c- 
Note that by taking the auxiliary component of X and 
inserting its VEV, (Fx) = /, one recovers the usual gaug¬ 
ino mass terms. The goldstino or sgoldstino interactions 
are obtained by taking the fermion or scalar component 
of X. 

We will from hereon focus on the interactions of the 
sgoldstino x = {(j) + ia)f^/2, where </> and a are the CP- 
even and CP-odd real scalar components. All the rele¬ 
vant vertices arising from eq. (2) can now be collected 
and included in the following sgoldstino Lagrangian [39] 


c _ 

®®" 2 v/ 2 / 


^ — ^gg + ^11 + + Fzz + Fww + > (3) 

where 

. + aW+ W, 

+ (4c) 

+ (4d) 


Lww = 


m 2 

717 ' 


I 1 




Fzz = 




2V2f 


= 


micl +17123^ 


9w 


Fz-i = 


2V2f 

{m2-mi)se„ce 


(4a) 

(4b) 


^^GG = 


2 ^ 2 / 


n^se„ce„ aF >^’', (4e) 

(-<()GG-biaGG)-bh.c., (4f) 


where sg„ = sui6w and ce„ = cos6w, with 9w be¬ 
ing the weak mixing angle, and the tilde denotes 
e.g. G“j, = (l/2)e^,ypcrG“^‘^. The interactions in eq. (4f) 
arise from the operator m^/(4/^)(XlX)^ in the Kahler 
potential, from which also the soft mass m^ = m^ for the 
CP-even and CP-odd sgoldstino scalars (j) and a arises. 
Notice that the sgoldstino couples purely to the trans¬ 
verse components of the W and Z bosons. A small cou¬ 
pling to the longitudinal components can arise through 
mixing with the Higgs, but for the region of parameter 
space that we consider such a mixing is negligible. 

From the sgoldstino Lagrangian (3) we can compute 
the partial decay widths for the sgoldstino scalar (p (the 


corresponding widths for a are obtained by simply re¬ 
placing (/)—>■ a since they are identical to those of (p), 
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(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

(5f) 


where the function k{x) = (1—4x2-|-6a;‘^)(l — is 

close to unity in the case where mw, fuz ^ m^ = 2 TeV. 

The interactions between the sgoldstino and the SM 
fermions arise from superpotential operators such as 
[Au/f)XQHuU‘^, which, upon taking the auxiliary com¬ 
ponent of X and inserting its VEV, also give rise to the 
usual A-terms. Since we are requiring all soft parameters 
to be smaller than y/J, the sgoldstino couplings will be 
suppressed at least by the ratio of the Higgs VEV over 
y/J, which makes the sgoldstino decays to SM fermions 
negligible with respect to the sgoldstino decays to vector 
bosons. 

The sgoldstino is produced at the LHC by gluon-gluon 
fusion with the leading order production cross section 
(summing the two equal contributions from the CP-even 
scalar p and the CP-odd scalar a, with m^ = ma) 


7 r2r {p -» gg) 
Asm^ 




( 6 ) 

where the partial width F (<() —>■ gg) is given by Eq. (5a), 
3 is the center of mass energy squared and fp/g (a:,Q^) 
are the parton distribution functions defined at the scale 
Q^. Since V^jm^ is below 10% for the sgoldstino in the 
relevant region of the parameter space, eq. (6), which 
assumes the narrow width approximation, is always reli¬ 
able. 


III. EXPLAINING THE DIBOSON EXCESS 

In this section, to assess the compatibility of a sgold¬ 
stino signal with the ATLAS diboson excess, we compare 
the number of signal events the sgoldstino gives rise to 
with the number of excess events reported by ATLAS. 
Figure 5 of ref. [1] shows the invariant mass distribution 
of the boson tagged jets for the tUZ, WW and ZZ selec¬ 
tion regions (SRs). These regions have large overlaps due 
to the overlap of the shapes of the single W and Z tagged 
jet mass distributions. We take into account this overlap 
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selection region 

final state 

WW 

WZ 

ZZ 

WW 

0.39 

0.37 

0.16 

WZ 

0.33 

0.44 

0.25 

ZZ 

0.27 

0.47 

0.37 


TABLE I. Efficiencies ^VfVj^VsV for ^ final state V/Vj to 
end up in the VsVi selection region. 


by computing the different efficiencies evfV^^VsV^ for a 
final state VfVj to end up in the VsV!. SR. The values of 
these efficiencies are given in Table I. 

We consider the window 1.75 — 2.25 TeV in the boson 
tagged dijet mass distribution, around the mass hypoth¬ 
esis, and compare the number of observed events with 
the number of events predicted by the SM. A reliable 
combination of the three SRs WZ, WW and ZZ would 
require detailed knowledge of the degree of correlation 
between these three channels that go beyond the effect 
of the efficiencies e that we take into account, such as the 
correlation of all the systematic uncertainties. Since we 
do not have this information at our disposal, we instead 
extract the signal from a single channel and then con¬ 
front it with the number of events observed in the other 
two channels, as well as with the other relevant analyses. 

Our model predicts the largest production rate in the 
WW channel and therefore, to also minimize the uncer¬ 
tainties coming from our extraction of the tagging effi¬ 
ciencies, we extract the sgoldstino signal yield from the 
WW SR. For an invariant mass of 2 TeV, we obtain from 
the ATLAS analysis Sww = 4 . 2 ^ 2 ’o excess events in the 
considered window, where the error band represents Pois- 
sonian central 68 % CL interval, which is what we refer 
to as Icr interval throughout the paper. The number of 
signal events produced by the sgoldstino in the various 
VsVg SRs is given by 


^v^vi = 


(TX BRwwx^ww X BRwiv-s-had X eww^v^v 


-I- crx BRzz X Azz X BR22_j.had X ezz^VsV' 




( 7 ) 

where BR^ corresponds to the sgoldstino decay branch¬ 
ing ratio into the ij final state, the factors Aww and 
Azz are the acceptances for the kinematic and topology 
selections and include the signal acceptance to the invari¬ 
ant mass cut in the window we consider, BRww->-had and 
BRzz-t-had are the hadronic branching ratios of the WW 
and ZZ channels, respectively, and C = 20.3 fb“^ is the 
integrated luminosity. 

The knowledge of the acceptance factors is a key in¬ 
gredient to estimate the number of events starting from 
a certain theoretical value of cr x BR. Unfortunately, 
the ATLAS analysis only reports the value of these ac¬ 
ceptances for a vector, and a spin-two signal hypothesis, 
the bulk graviton, that decays into longitudinally polar¬ 
ized vector bosons. This information does not allow us 
to extract the acceptance for a scalar particle decaying 


to transverse vector bosons, as is the case of the sgold¬ 
stino. Comparing the ATLAS acceptances with the ones 
that CMS reports in ref. [2], which is the counterpart of 
ref. [ 1 ], where also a spin-two particle decaying to trans¬ 
verse vector bosons is considered, the RS graviton, we ex¬ 
pect the acceptances for a resonance decaying into trans¬ 
verse gauge bosons to be about 50% smaller than the 
acceptances for a resonance decaying into longitudinal 
vectors, which is the case reported by ATLAS. Such a re¬ 
duction in the acceptances would require a larger a x BR 
to explain the ATLAS excess, thereby selecting a region 
of the parameter space with lower values of the relevant 
parameters, namely m 2 , m 3 and ^/J. However, since we 
cannot reliably estimate this number we will instead use 
the acceptances that ATLAS reports for the bulk gravi¬ 
ton also for the sgoldstino. This also allows us, in a more 
consistent way, to compare with other analyses where 
the same spin-two signal hypothesis is considered. We 
stress that, in the case in which a higher significance of 
this excess is observed in Run II, it is of primary impor¬ 
tance that the experimental collaborations provide the 
acceptances for all the relevant spin hypotheses and po¬ 
larizations of the final state vector bosons. 

The acceptance factors in eq. (7) are estimated start¬ 
ing from the total selection efficiencies reported in Fig¬ 
ure 2 (b) of ref. [ 1 ] , divided by the aforementioned boson 
tagging efficiencies in the respective SRs. This number 
is then multiplied by the efficiency corresponding to the 
invariant mass cut in the window we consider, estimated 
from the signal shape reported by ATLAS in the WW 
SR for the spin-two resonance in Figure 5 (b), which is 
the one with more available statistics and which is about 
0.87. The resulting two acceptances in eq. (7) are, as ex¬ 
pected, almost identical, Aww = 0-22 and Azz = 0 . 21 . 

The parameters that most strongly affect the diboson 
channels relevant for the excess are m 2 , m 3 and , while 
in the case of the 77 and Z 7 channels, there is also some 
dependence on mi. The ATLAS search in the 77 channel 
of ref. [50] place a 95% CL upper limit on cr x BR-^.^ at 
around 0.3 fb for a mass of 2 TeV. We stress that also 
this result is obtained assuming a spin-two resonance and 
can not be straightforwardly used to constrain a scalar. 
However, since we do not expect huge changes in the ef¬ 
ficiencies, this gives us an estimate of the current bound 
on a heavy scalar decaying to 77 . As can be seen from 
eq. (5d), this bound can be completely evaded by choos¬ 
ing mi « —m 2 tan^ 9w since, in this case, BR.y.y « 0 . 
However, since we find that there is a wide range of mi 
that satisfies the 77 constraint, without affecting the di¬ 
boson channels, we choose not to fix any particular re¬ 
lation to m 2 , but instead we set it to a reference value, 
mi = 100 GeV^, for which the constraint from the 77 
search is satisfied in the entire range or m 2 , m 3 and i/J 


^ Notice that when mi < ra^jl the sgoldstino can decay into 
two neutralinos with a coupling that is generally model depen¬ 
dent. This would amount in a shift to the total width that would 
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FIG. 1. The four regions enclosed by the solid lines, corre¬ 
sponding to four different values of m 2 , show the points in 
the plane (m 3 , v7) for which a 2TeV sgoldstino can explain 
the ATLAS diboson excess [1] within la. The (dark blue) 
regions enclosed by the dashed lines correspond to the subset 
of these point that satisfies the 95% CL limit placed by the 
semi-leptonic diboson search by CMS [51]. 

that we consider. 

The regions enclosed by the solid grey curves in Fig¬ 
ure 1 represent the regions in the (m 3 , ■\/f) plane where 
the 2 TeV sgoldstino signal reproduces the number of ex¬ 
cess events observed by ATLAS in the WW SR in the 
invariant mass window we consider, within the la band, 
namely Sww = 4.2l2'o- different regions cor¬ 

respond to four representative values of m 2 , namely 1 , 
2.5, 4 and 5 TeV. Of course, different values of m 2 in¬ 
terpolate between these regions covering a large part of 
the m 3 < vT plane. In Figure 1 we have also imposed 
the constraint m 2 , m 3 < v 7 to ensure a valid expansion 
in eq. (3) in terms of effective operators. The current 
experimental limits on m 3 and \/~f depend on the masses 
of the other superpartners on which the fit to the AT¬ 
LAS excess imposes no constraints. In Figure 1, we take 
a conservative approach and require both m 3 and v 7 
to be above 1.5 TeV, corresponding to the most strin¬ 
gent current bound on both the gluino mass [52] and the 
SUSY breaking scale [38]. 

In the four regions shown in Figure I, the predicted 
values for the cross sections in the other two diboson 


slightly shift the region of the parameter space where the model 
reproduces the ATLAS excess. This decay could also give rise 
to a final state with two photons and missing transverse energy 
that is potentially interesting at the LHC. 


channels span the ranges a x BRwz G [2.3, 7.6] fb and 
the cr X BRzz G [0.7,2.2] fb. For these cross section in¬ 
tervals, by using eq. (7), we compute the Icr intervals for 
the number of excess events predicted in the WZ and ZZ 
SRs, and we obtain Swz G [2.4,8.0] and Szz G [1.3,4.2]. 
Comparing these numbers with the corresponding val¬ 
ues extracted from the ATLAS analysis, namely 7 . 0 l 2'6 
and 6 .4:^2A’ respectively, we see that the WZ SR is well 
within the statistical Icr band, while ZZ shows a slight 
tension. However, this tension is removed once one in¬ 
cludes systematic uncertainties, which are at the level of 
50% for the signal [1]. Nevertheless, it is worth asking if 
we can directly understand this tension from the ATLAS 
analysis. The ZZ SR in the ATLAS analysis is the one 
suffering from less statistics and is therefore expected to 
have the largest uncertainty. One important feature that 
one notices looking carefully at the ATLAS analysis is 
that the shape of the fitted background in the WZ and 
WW channels agree rather well, while the one in the ZZ 
channel has a different shape at high masses, as it falls 
much faster. Since we do not expect substantially dif¬ 
ferent shapes for these SM diboson backgrounds at such 
high invariant masses, this provides an estimate for the 
error of the fit to the background. In order to assess the 
possible origin of the tension that we find between the sig¬ 
nal in the ZZ SR and in the other two SRs, we compute 
the number of excess events in the ZZ SR that would be 
obtained if we instead use the WZ ox the WW shape, 
with the ZZ normalization, as ZZ background distribu¬ 
tion. With the WZ ox the WW background distribution 
shapes, we obtain only 3.0 and 1.7 excess events, respec¬ 
tively, in the window we consider, to be compared with 
the 6.4 obtained from the ATLAS fit to the background 
distribution in the ZZ SR. This shows that the tension in 
the ZZ channel could be due to an inaccurate prediction 
of the shape of the background in the ZZ SR, caused by 
the lack of statistics. 

IV. OTHER ANALYSES 

Now that we have extracted the interesting region of 
the parameter space of the sgoldstino that allows us 
to reproduce the excess of events observed by ATLAS, 
we confront our signal hypothesis with the other rele¬ 
vant searches. The first search to compare with is re¬ 
ported in ref. [2] and is the CMS analogous of the AT¬ 
LAS fully hadronic search [1]. This sets the limits (for 
the same spin and polarization hypothesis that is consid¬ 
ered by ATLAS and that we used to extract the signal) 
a X BRww < llfb and a x BRzz < 10 fb for a 2 TeV 
mass hypothesis, which lie above our la bands for the 
corresponding quantities and hence do not set any fur¬ 
ther constraint on the allowed parameter space. 

The analysis that sets the strongest constraint on the 
a X BR for a 2 TeV resonance decaying to gauge bosons 
is a CMS search in the semi-leptonic channel with either 
1 lepton (WW channel) or 2 leptons (ZZ channel) [51]. 
There is no overlap in this case due to the selection with 
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different numbers of leptons. We can therefore directly 
compare our predicted a x BR with the limits this search 
places in the respective channels, tr x BRiyvv < 3 fb and 
a X BR^z < 8 fb at 2 TeV. While the ZZ bound does 
not constrain our parameter space, we get a constraint 
from the WW channel, which reduces the allowed pa¬ 
rameter space in the plane in Figure 1 to the 

(dark blue) regions enclosed by the dashed lines. The 
(light blue) regions that remain outside the dashed con¬ 
tours are excluded at 95% CL by the CMS semi-leptonic 
analysis [51] in the WW channel. 

Let us finally comment on other possible interesting 
channels. If the ATLAS diboson excess is caused by the 
scalar sgoldstino, no signal is expected in the ZH and 
WH channels. Hence, if statistically significant excesses 
are found in Run II in these channels, it would point 
toward other new physics scenarios. Instead, the most 
relevant other channels for the sgoldstino signal hypoth¬ 
esis are the 77 and Zj channels. As was discussed in the 
previous section, the bound from existing 77 searches 
can always be satisfied by choosing mi to be within a 
rather wide range around the value mi « —m 2 tan^ 
for which BR.^.^ vanishes, as can be seen from eq. (5d). 
Clearly, the allowed range of mi is wider for larger values 
of m 2 . 

One way to place a constraint on mi would be to 
search for a resonance in the Z 7 channel at 2 TeV, and 
use the relation between mi and m 2 in eq. (5e). The 
only search in the Z 7 channel that we are aware of is the 
ATLAS analysis in ref. [53] which only extends to invari¬ 
ant masses up to 1.6 TeV. Therefore we do not get any 
constraint from this search. However, it is interesting to 
note that the exclusion at 1.6 TeV is u x BR^.^ < 0.17 fb 
for a scalar. We find that, in most of the the parameter 
space that explains the ATLAS diboson excess, and for a 
large range of values of mi within the region allowed by 
77 searches, the a x BRz-y we get for the 2 TeV sgold¬ 
stino is larger than 0.17fb. This suggests that once the 
Z 7 search is extended to include 2 TeV invariant masses, 
this channel will be sensitive to the sgoldstino signal and 
could quickly lead to a discovery or exclude most of the 
parameter space presently allowed. 

V. CONCLUSIONS 

In this paper we have provided an explanation of the 
recently reported ATLAS diboson excess in terms of a 


2 TeV sgoldstino scalar, which is present in a class of su¬ 
persymmetric models in which the supersymmetry break¬ 
ing scale is in the few TeV range. Fitting this excess se¬ 
lects particular ranges and relations among the gaugino 
masses, while imposing no constraints on the other su¬ 
perpartner masses. In terms of other resonance searches, 
while no signal is expected in the ZH and WH channels, 
we expect the most sensitive channel to be Z 7 . 

The sgoldstino production cross-section, which origi¬ 
nates from gluon-gluon fusion, is expected to increase by 
a factor of about 19 when going from = 8 TeV to 
13 TeV for a mass of 2 TeV. This should be contrasted 
with, for example, the factor of about 7 increase of the 
production cross-section that is expected for a qq reso¬ 
nance of the same mass. Hence, with the sgoldstino signal 
hypothesis, also taking into account that the background 
is mainly due to qq, one expects the diboson excess to 
grow significantly faster with the incoming 13 TeV data, 
with respect to, for instance, a heavy vector signal hy¬ 
pothesis. The different scaling of the signal cross sections 
with the collider energy could help, in case of discovery, 
to understand the nature of the new resonance. 

Let us end by encouraging the ATLAS and CMS col¬ 
laborations to provide the efficiencies for all the relevant 
spin hypotheses of the resonance and the polarizations 
of the vector bosons in the final state. The fact that 
the ATLAS analysis [1] only provides the efficiencies for 
spin-one and -two resonances decaying to longitudinally 
polarized vector bosons introduces a large uncertainty in 
our interpretation of the excess. However, we expect this 
to only amount to a rescaling and possibly a shift of the 
relevant parameter space region towards slightly lower 
values of m 2 , m 3 and v 7 - 
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